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Abstract 

The paper investigates possibility of equilibrium solid-body rotation of a vortex bundle diverging 
at some height from a cylinder axis and terminating on a lateral wall of a container. Such a bundle 
arises when vorticity expands up from a container bottom eventually filling the whole container. 
The analysis starts from a single vortex, then goes to a vortex sheet, and finally addresses a multi- 
layered crystal vortex bundle. The equilibrium solid-body rotation of the vortex bundle requires 
that the thermodynamic potentials in the vortex-filled and in the vortex-free parts of the container 
are equal providing the absence of a force on the vortex front separating the two parts. The 
paper considers also a weakly non-equilibrium state when the bundle and the container rotate with 
different angular velocities and the vortex front propagates with the velocity determined by friction 
between vortices and the container or the normal liquid moving together with the container. 

PACS numbers: 67.30.hb,47.15.ki, 67.30.he 
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I. INTRODUCTION 



Transient processes of establishing of stable vorticity in rotating superfluids have always 
been in the focus of attention of theorists and experimentalists studying superfluid vortex 
dynamics. An important example of such a process is penetration of a vortex bundle into an 
originally vortex-free rotating container filled with a superfluid. This process was thoroughly 
investigated in superfluid ^He-B theoretically and experimentally-"-. Vorticity is generated 
at the container bottom and propagates upward along the cylindric container axis in the form 
of a vortex bundle flaring to lateral container walls. The flaring part of the vortex bundle was 
called vortex front. Below the front the vortex bundle is vertical but twisted. The twist is 
connected with the flux of the angular momentum along the bundle, which must dissipate due 
to either mutual friction in the bulk or friction at rough wall surface. A great attention was 
directed to transition from laminar to turbulent vortex front propagation, especially at low 
temperature where disappearance of mutual friction facilitates the transition to turbulence. 

All studies of the flaring vortex bundle, analytical, numerical, and experimental, were 
performed in the presence of dissipation, without which front propagation is impossible 
since it is accompanied by change of the total energy and the angular momentum. The goal 
of the present work was to check whether a stable solution for a vortex bundle terminating 
on a lateral wall may exist as an equilibrium state without dissipation and propagation 
along the rotation axis. In this equilibrium state the whole bundle together with its vortex 
front rotates without twisting as a solid with constant angular velocity. If the container 
rotates with the same angular velocity neither dissipation nor propagation of the vortex 
front along the rotation axis is possible. Our paper considers conditions for existence of 
such "eigenrotation" and analyses the effect of weak friction, which leads to propagation of 
the vortex front if the container and the bundle rotate with different angular velocities. 

The paper starts from the predecessor of the vortex bundle: a single vortex filament, 
located on the container axis in the lower part of the container, at its higher part continuously 
goes away from the axis and eventually terminates on the lateral wall (Sec. HI])- Section HIT] 
analyses a vortex bundle in which vortices form a single axisymmetric layer (vortex sheet). 
Section IIVI addresses the case of two coaxial non-interacting sheets and demonstrates that 
it is impossible to find the equilibrium solid-body rotation for such a two-layer bundle. 
This led to conclusion that equilibrium solid-body rotation of a multi-layered vortex bundle 
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FIG. 1. (Color online) Vortex line attached to a thin wire. The profile of the vortex line can be 
examined on two different length scales: (a) At large scales r ^ the local induction approxi- 
mation becomes applicable, in which the wire is treated as an enhanced vortex core. The vortex 
line meets a wire at the finite contact angle 9. (b) On small scales r ^ the vortex line goes 
smoothly over into the image vortex (dashed continuation of the vortex line) , with the wire surface 
perpendicular at the connection point. 

requires effective interaction between layers. This could be the same interaction, which leads 
to formation of a vortex crystal. The condition of equilibrium rotation of a stiff solid-body 
vortex bundle diverging to a lateral wall at some height is analyzed in Sec. |Vl Section I VI I 
discusses weakly non-equilibrium vortex bundle, when friction of the rotating vortex bundle 
makes possible propagation of the vortex front along the rotation axis. Concluding remarks 
are presented in the last Sec. IVIIi 

II. SINGLE VORTEX 

Let us start from the most elementary case when the bundle reduces to a single vortex, 
which terminates on a lateral wall of a container. A freely precessing vortex is in fact a 
particular case of the geometry, which has already been carefully investigated theoretically 
and experimentally in connection with the study of precession of single vortex trapped on 
a wire coaxial with a cylindric container filled by superfluid ^He-B^"®. Geometry of the 
experiment^ is shown in Fig. [TJ^a). The z-axis is the axis of the cell of the radius R and 
of the wire of the radius r^. The vortex filament with the circulation n of the superfluid 
velocity being trapped on the wire at z < Zq, peels from the wire at z = Zq stretching to 
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the container lateral wall. Just the "unzipped" (free) part of the vortex filament aX z > zq 
participates in the precession. 

If dissipation is absent, the vortex-line shape may be found by minimization of the energy 
in the coordinate frame rotating with the angular velocity Vt of the vortex precession. This 
energy corresponds to the Gibbs thermodynamic potential 



G = E- QM, 
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is the energy in the local induction approximation, 
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is the angular momentum of the liquid around the axis z, is the vortex core radius, and 
the shape of the vortex line in the cylindric coordinates is given by the function z{r). It is 
supposed that the line is inside the axial plane, so the azimuthal angle does not vary along 
the line. The function z{r) is determined by variation of the thermodynamic potential with 
respect to the vertical coordinate z of the vortex line. The corresponding Euler-Lagrange 
equation is 

e{r-) dz/dr 
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Att r^ 

is the energy per unit length of the vortex line, which determines the line-tension force. 
Equation (jl]) was derived from minimization of the thermodynamic potential, but at the 
same time it is the equation of balance of forces on the precessing vortex line: the Magnus 
force on the left-hand side, which is determined by the linear velocity Qr of the precession, 
is balanced by the line tension force oc e. One should look for a solution of this equation 
with the boundary conditions on the lateral wall and the wire. The first one is the condition 
of transversality (the vortex line normally ends at the lateral wall of the cylinder): 
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The second boundary condition is imposed in the "unzipping" point z = zq'- 
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where 9 is the angle between the z axis and the vortex hne at z = and = {ps^^/ 47r) In (R/ r 
is the hne tension of the vortex filament attached to the wire. The condition directly fol- 
lows from minimization of the Gibbs potential with respect to variation of the unzipping 
coordinate It ensures the balance of line-tension forces along the z axis for the element 
of the vortex filament in the point z = Zq (see Fig. [1]): The unzipping point is at rest if 
Ew = £ COS 9. In our approach the trapped segment of the vortex line is treated as a vortex 
line with a larger core radius r^,. So the approach is valid only for rather thin wires with 
radius much less than the radius R of the container. 

Integrating Eq. (jl]) over r with the transversality condition Eq. 
integral 



one obtains the first 
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At the unzipping point this equation together with boundary condition Eq. ((Tj) yields that 



g = e — flniz = 0, 



(9) 



where g, e = e^, and = psnfl{r'^ — r1/)/2 are the Gibbs potential, the energy and the 
angular momentum per unit length of the vortex line below the unzipping point. Emergence 
of the condition of zero Gibbs potential below the unzipping point from the minimization 
of the total Gibbs potential is quite natural. In the vortex-free region the Gibbs potential 
density vanishes, and if the Gibbs potential density below the unzipping point is nonzero, 
then there is a force driving the curved piece of the vortex line (an analog of the vortex 
front) upward or downward. The condition ([9]) determines the precession angular velocity 
of the curved vortex line stretched between the wire and the cell lateral wall^: 
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Note that this expression is exact and does not depend on vortex line shape of the curved 
segment or on using the local induction approximation^"-. This is a direct consequence of 
the rigorous canonic relation 
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Varying the position of the unzipping point the variations of the energy and momentum 
are dE = Swdz^ and dMz = rrizdzo, so OE/dMz = Sw/mz, and Eq. f lTTj) yields expression 
Eq. fllOp for the precession angular frequency. 

This analysis can be applied to the case of a free axial vortex ending at the wall. The 
case is the limit of an extremely thin wire when the wire radius must be replaced by the 
vortex core radius in all expression. Formally this means that the unzipping point goes to 
z —7- — oo, and the curved vortex line smoothly approaches to the vertical axis of rotation in 
accordance with the boundary condition Eq. ([7]), which now tells that 
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The shape of the free vortex can be analytically obtained after the second integration of 
Eq. ([8]) taking into account the expression Eq. (fTOj) for the angular velocity Q: 

zir) ^ V2W^ -R-^ln ^^^^f^\ (13) 

The expression demonstrates that the vortex line exponentially approaches to the axis r = 0: 

While the frequency of the single vortex precession straightforwardly follows from com- 
monly accepted thermodynamic arguments, the shape of the vortex line was a matter of 
dispute resulting from disagreement on a proper usage of the local induction approximation 
for the precessing partially trapped vortex line. One may find a detailed discussion of the 
issue m Refs. HandS In particular, the debate was about a proper boundary condition 
at the unzipping point. Instead of the boundary condition Eq. ((Tj) based on the balance 
of forces directly following from the variational principle Schwarz- used the condition that 
the vortex line is normal to the wire. At the very surface of the wire the latter boundary 
condition is definitely correct. But at small scales of the order of the wire radius there are 
forces, which led to fast deviation from the normal direction. These forces were neglected 
by Schwarz and all others addressing this problem. This is a legitimate approximation when 
one looks for the vortex shape at large scales of the order of the container radius R but 
only if one uses the boundary condition ([7]) based on the balance of line-tension forces. In 
reality this means that the boundary condition is imposed not exactly at the radius of 
the wire but on the distance larger than r^, which at the same time is still much smaller 
than i?, as illustrated in Fig. [T|d. It is worthwhile to note that the analysis of the shape of a 



free precessing vortex gives one more justification of the force-balance boundary condition. 
The latter provides a natural transition from the vortex partially trapped by the wire to a 
free vortex smoothly changing its direction from vertical to horizontal. On the other hand, 
Schwarz's boundary condition becomes senseless for the free vortex since it requires that the 
vortex meets the axis normally. 



III. SINGLE VORTEX SHEET 



Our next step is to analyze a bundle of vortices but still of simple geometry: A^i vortices 
form a cylindric vertical sheet of radius ri which at some height diverges to lateral wall 
forming a whorl (see Fig. 2). The single-sheet whorl is a simulation of a more complicated 
vortex front. As in the case of a single vortex, any vortex line in the sheet is given by the 
function z{r), which is independent on (p and is determined from the variational principle 
for the the total Gibbs potential G = E — QMz- The total angular momentum is 
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where Vs{r) = Nin/2nr is the azimuthal superfiuid velocity induced by the vortex sheet. 
The total kinetic energy E = Eg + consists of the energy 
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of the velocity field induced by the vortex sheet and the energy of individual vortex lines 
given by [compare with the expression Eq. ([2]) for a single vortex]: 



An Jn Tr 



\ 



1 + 



dz{r) 
dr 



1 2 



dr. 



(16) 



where b = 2nr/Ni the r-dependent intervortex spacing. 

Variation of the Gibbs potential of the vortex sheet yields the Euler-Lagrange equation: 

\2 1 AT. . ..2 ^ , X dz/dr 



p^ 



^ PsiNinY 1 NiPsk' d (2nr\ 

:NiKUr = — —In — — 

An r An dr KJMirc/ h 



(17) 



^l + {dz/drY' 

The first integral of this equation for the boundary condition Eq. (E]) at the container lateral 
wall is 
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FIG. 2. (Color online) Single vortex sheet diverging to lateral walls via forming a whorl 



The boundary condition 



dz 
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r=ri 



provides a transition from the whorl to the vertical stem of the vortex sheet. Using it in 
Eq. (fT8|) gives the condition that the Gibbs potential per unit length gi = e — QrUz {e and 
rriz are the energy and the angular momentum per unit length of the stem) vanishes at the 
stem: 
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As in the case of a single vortex, this condition is necessary for the absence of the force 
driving the front along the vertical axis. 

But our minimization of the Gibbs potential requires an additional step: minimization g 
with respect to ri at fixed number of vortices. The minimization yields that 
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The relation means that the whole sheet rotates as a solid body with linear velocity Qri = 
{Ni — l)K,/ATTri. Though our analysis was based on the local induction approximation, for the 
vertical part (stem) of the vortex sheet the obtained expression is exact and is easily derived 
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from the Bio-Savart law for A^i equidistant vortices on the circumference of the radius ri. 
The modulus of the velocity induced by one vortex on another, Vmd = k/2tt6, is determined 
by the distance 6 = 2ri sin(a/2) between them, where a is the angle between radii directed 
from the center to the both vortices. The azimuthal component f sin(a/2) = n/ATrri of 
the velocity does not depend on distance between two vortices, and therefore any vortex in 
the sheet moves along the circumference with the azimuthal velocity v sheet = {N — Vjn/Airri 
induced by rest — 1 vortices. The radial component of any vortex vanishes for equidistant 
location of vortices by symmetry, li Ni ^ 1 the velocity v sheet = {vin + Vout) is an 
average of the velocities fj„ = and Vout = NiK/2Tiri on both sides of the sheet (inside and 
outside). This law of motion for vortex sheets is well known in classical hydrodynamics. 

In Fig. |3]we depicted the solution of the Euler-Lagrange equation ( !T7|) with the boundary 
conditions z{R) = and dz/drln = for A^i = 10 and QR^/k = 4.52828. The whorl 
smoothly diverges from the vertical stem with the radius determined from Eq. (12T]) and 
terminates on the lateral wall. 



IV. TWO VORTEX SHEETS 

One might think that a more realistic vortex bundle could be modeled as an ensemble 
of coaxial vortex sheets with their stems rotating together with the same angular velocity 
as a solid body. In order to check this option we considered two coaxial vortex sheets with 
numbers of vortices A^^i and N2 diverging to lateral walls via two whorls (Fig. H]). The 
analysis of the Gibbs potential of the sheets below the whorls has shown that the condition 
for solid body rotation of the two sheets with the same angular velocity cannot be realized. 
The density of the Gibbs potential for two vertical sheet is g = gi + g2 where the Gibbs 
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FIG. 3. The profile of the vortex sheet calculated from the Euler-Lagrange equation ()17p for 
Ni = 10 and CIR'^/k = 4.52828. 
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FIG. 4. (Color online) Hypothetical picture of two vortex sheets terminating on the lateral wall. 
But no state of solid-body rotation of two sheets with the same angular velocity was found. 

potential gi per unit length of the inner sheet is given by Eq. (120|) and the Gibbs potential 
g2 per unit length of the outer sheet is 

P..'N,i2N, + N,) /fix iW ,^ (2^^ _ _ .^ 



Minimization with respect to r2 yields the relation 

{N2 + 2Ni - 1)k 



n 



(23) 



The equilibrium condition (absence of the driving force on the whorl) is satisfied only if (71 = 
and g2 = 0. One may find the states with vanishing gi and g2 only if the angular velocities 
Q in Eqs. ( 12T]) and ( l23l) are different. The state with the solid body rotation of two sheets 
together cannot be found. However this outcome results from shortcoming of our model and 
does not mean that the solid-body "stem + whorl" structure for multi-layer vortex bundles 
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is impossible. In our model the vortex sheet induces a fully axisymmetric velocity field 
outside the sheet and no field inside. Though discrete vortex structure of the sheet is taken 
into account via a logarithm contribution in the energy, the effect of individual vortices on 
the velocity field is neglected. Such an approximation is valid only if the distance between 
sheets exceeds the intervortex distance in any sheet. In a tightly packed vortex bundle 
the approximation fails, and interaction between neighboring layers glues them effectively 
and does not allow rotation with different velocities. Then the condition for an equilibrium 
vortex front becomes much less severe: it is not necessary for the Gibbs potential of any 
sheet to vanish, it is sufficient that the total Gibbs potential vanishes. The latter condition 
will be considered in the next section. 



V. EQUILIBRIUM ROTATION OF THE SOLID-BODY VORTEX BUNDLE TER- 
MINATING ON THE LATERAL WALL 

From the analysis of a single and two vortex sheets it is evident that the most restrictive 
condition for equilibrium solid-body rotation of a vortex bundle terminating on the lateral 
wall is zero Gibbs potential much below the front. If this condition is satisfied the solution 
of the differential Euler-Lagrange equation for the front shape is straightforward, though 
technically complicated: for a multilayered vortex whorl one should solve a partial differential 
equation in the space of two coordinates r and z. This section addresses only the condition 
of solvability of this equation: absence of a force on the vortex front. 

We look for the equilibrium state of fixed number of vortices forming the solid-body 
vortex bundle terminating on the lateral wall and rotating with the angular velocity Q (Fig. 
[5]). The azimuthal velocity field in the stem of the bundle is 
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where Rq is the radius of the bundle stem. The number of vortices and the bundle radius 
are connected with the relation = 27tQRq/k. The Gibbs potential for the unit length of 
the bundle stem in the coordinate frame rotating with the angular velocity Qq is 

g = e- n,m, = np^n'R', fin ^ + 7) + ^^^^ In ^ - fiovrp.fii?^ ( R^ - ^] . (25) 
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FIG. 5. (Color online) Solid-body vortex bundle terminating on the lateral wall. 

In the equilibrium Q must coincide with Qq and using the relation between Q and Rq this 
equation reduces to 

,2 
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The condition g = yields the relation 



(27) 



If the bundle nearly fills the whole container [Rq ^ R, i.e., the number N of vortices is close 
to the equilibrium value) 

In - f« — + In y/N. (28) 

So the condition for equilibrium rotation of the bundle terminating on the lateral wall can 
be satisfied but not for a very large number of vortices. If the filling factor Rq/R is small, 
unrealistically high values of ln(i?/ Vc) are required for equilibrium rotation of a large number 
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of vortices. Equilibrium rotation is impossible in the classical limit of continuous vorticity 
neglecting the vortex line tension oc «;ln(r^/rc). 

VI. NON-EQUILIBRIUM ROTATION AND FRICTION 

The equilibrium finite vortex bundle may rotate freely only if a container rotates with the 
same angular velocity or there is no friction force between moving vortices and the container 
and the normal liquid rotating with it. In the experiment this condition is not satisfied, 
and dynamics and structure of the bundle require a more complicated analysis in general. 
But if friction is rather weak as expected in the low-temperature limit, one may consider 
the friction effect assuming that the structure of the vortex front (whorl) is not affected 
seriously. 

For the state close to the equilibrium one may use the thermodynamic approach. If the 
vortex front moves with the velocity Vf variation of the energy of the vortex bundle rotating 
with the angular velocity during a short time interval dt is (the energy of the front is 
neglected compared to that of the bundle stem) 

dEi, = e{Q)vfdt. (29) 

Displacement of the vortex front is accompanied by variation of the bundle angular momen- 
tum dMz — rrizVfdt. The total angular momentum is conserved, and the angular momentum 
dMz must be transferred from the container via mutual friction with the normal liquid mov- 
ing rigidly with the container or via surface friction of vortex ends at a rough wall. This 
leads to energy variation of the container rotating with the angular velocity Qq'- 

dEc = nodM^. (30) 

The total decrease of the energy due to front motion is 

dE = dEc - dEb = -[e(Q) - nQm^{n)]vfdt. (31) 

This supports the condition of the equilibrium at fl — flo used above: the Gibbs potential 
density g{ft) — e(Q) — ^rrizift) must vanish. The quantity Ff — e(Q) — ilouizi^) can be 
interpreted as a force driving the vortex front in a non-equilibrium state. Since the total 
energy cannot change the released energy dE is compensated by dissipation with the rate 
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Q = dE /dt = FfVf. Eltsov et al- estimated the dissipation rate Q assuming that the bundle 
with Rq ^ R rotates with the same velocity as the container {Q = Qq) and neglecting the 
quantum line-tension contribution oc K,ln{ry/rc) to the energy. This yielded the force Ff 
equal to the bundle kinetic energy ixpsVl^R'^ / ^ per unit length. Here we consider the state 
close to the equilibrium solid-body rotation at which the force Ff = e{Q)~Qmz{^l) vanishes. 
So the force is 

Ff = [il - llo}m^{n) = . (32) 

The vortex front can move only if there is exchange of the angular momentum between 
the bundle and the container. We assume the phenomenological linear relation between the 
friction torque t = '-/{Q — Qq) on a. vortex and the relative angular velocity f2 — i^o- Then 
the balance equation for the total angular momentum of the liquid is 

dM, 



m 



.Vf = Nt = ^N{Q-Qo). (33) 



dt 

For the vortex bundle with Rq ^ R the angular momentum per unit length is 



PsNkR^ /4, and the front velocity 



does not depend on the number of vortices in the bundle at given f2 — r2o. 



VII. DISCUSSION AND CONCLUSIONS 



The paper addressed the question whether the vortex bundle terminating on the lateral 
wall can rotate as a solid body and what is the angular velocity of such "eigenrotation" . 
There are two conditions for existence of eigenrotation: (i) The whole vortex bundle includ- 
ing its stem and the whorl is in an equilibrium solid-body rotation with the same speed; 
(ii) The Gibbs potential per unit length of the stem is equal to the zero Gibbs potential 
of the vortex-free state above the bundle. The latter condition eventually determines the 
angular velocity of the eigenrotation and was used in the past for estimation of the rotation 
of the vortex front by Eltsov et al}-^ but without paying attention to the first condition 
of equilibrium solid-body rotation of the whole bundle. Indeed, using Eq. (!25ll for linear 
density g of the Gibbs potential with the filling factor Rq/R ~ 1 and the logarithm term 
cx ln(r„/rc) neglected the condition g = yields that Qq = Q/2. So in the coordinate frame 
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rotating with this angular velocity the Gibbs potentials above and below the vortex front 
are equal, and on the basis of it Eltsov et al. concluded that the front must rotate with the 
speed, which is half of the rotation speed of the container. This rotation has nothing to do 
with the equilibrium rotation analyzed in the present work: The vortex front rotates twice 
slower than the vortex bundle stem, which must leads to steady twisting of the vortex bun- 
dle. Application of thermodynamic-balance arguments to so strongly non-equilibrium state 
requires justification, and a more rigorous dynamical approach of the problem is wanted. 
The analysis of the equilibrium and weakly non-equilibrium vortex bundle presented in the 
present paper can be considered as the first step in this direction. 
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